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Context. Optical spectra of classical T Tauri stars (cTTS) are rich in emission lines of low-excitation species that 
are composed of narrow and broad components, which indicates the existence of two emitting regions with different 
kinematics, densities, and temperatures. The photospheric spectrum is often veiled by an excess continuous emission. 
This veiling is usually attributed to radiation from a heated region beneath the accretion shock. The broad emission 
lines of Hi, Hen, Can, Fen, and other species are thought to form in a larger volume of gas. 

Aims. The aim of this research is to clarify the nature of the veiling, and whether the narrow chromospheric lines of 
Fe I and other metals represent a standard chromosphere of a late-type star, or are induced by mass accretion. 
Methods. We carried out high-resolution spectroscopy of selected cTTS with a special focus on DR Tauri and followed 
variations of chromospheric features, such as narrow Fel emission lines, and accretion signatures such as the veiling 
continuum and the He ll line emission. 

Results. We found that the amount of veiling in DR Tau varies from practically nothing to factors more than 10 times 
the stellar continuum intensity, and that the veiling is caused by both a non-photospheric continuum and chromospheric 
line emission filling in the photospheric absorption lines. The latter causes differential veiling because stronger lines are 
more veiled. We developed methods to separate the two sources of veiling. Several veiled T Tauri stars show a common 
effect: the radial velocities of photospheric and chromospheric lines vary in anti-phase. This is caused by an area with 
enhanced chromospheric emission, which is offset from the pole of rotation and is associated with the hot spot formed 
at the footprint of the magnetic funnel of mass accretion. 

Conclusions. The enhanced chromospheric emission in cTTS is linked not only to solar-like magnetic activity, but is 
powered to a greater extent by the accreting gas. We suggest that the area of enhanced chromospheric emission is 
induced by mass accretion, which modifies the local structure of stellar atmosphere in an area that is more extended 
than the hot accretion spot. The narrow emission lines from this extended area are responsible for the extra component 
in the veiling through line-filling of photospheric absorption lines. 
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1. Introduction 

The T Tauri stars are pre-main sequence stars of low mass. 
Spectroscopically they are characterized by emission lines 
superimposed on a photospheric spectrum of late type. In 
the classical T Tauri stars (cTTS) the emission line spec- 
trum is very rich and related to magnetospheric accretion 
from dusty disks. The T Tauri stars rotate faster than most 
late-type main-sequence stars, and the enhanced emission 
in the weak-line stars, those not exhibiting evidence of ac- 
cretion, is adequately explained in terms of enhanced solar- 
like magnetic activity. The same dynamo effect should be 
present in the cTTS as well, but in addition enhanced line 
emission is generated in the accreting gas, which flows from 
the circumstellar disk to the star, where it eventually heats 
up and releases energy in strong shocks at the stellar sur- 
face. Accordingly, the emission lines are much broader in 
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classical T Tauri stars than expected from chromospheric 
emission alone. For reviews of properties and models of 
cTTS see Petrov (|2003l) and Bouvier et al. p007|) . 

In cTTS the photospheric absorption lines are weakened 
by what is commonly attributed to an overlying continu- 
ous excess emission. The energy distribution of this veiling 
has been matched with black body radiation at a tempera- 
ture of ~ 10"* K, or alternatively with Paschen continuous 
emission (Basri & Batalha 1990; Hartigan et al. 1991, 19951 
Va lenti e t aL [T993l Calvet & GuUbring IT9981 Gullbring et 
al. I1998P . The veiling is attributed to radiation from "hot 
spots" , small areas (sometimes rings) of heated gas at the 
footprints of the accretion flows. T Tauri stars are pho- 
tometrically variable, both from rotational modulation of 
surface features and stochastic events (e.g., flares, changes 
in mass accretion rates) . The changes in line intensities and 
profiles in many cTTS are commonly attributed to varia- 
tions in the accretion flow. The veiling levels can vary dra- 
matically, and in the most active cTTS veiling may change 
from a tiny fraction of the intensity of the photospheric 
continuum to factors of several times this level. The degree 
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of veiling is commonly used to estimate accretion rates in 
CTTS. 

Intuitively, one would expect that when the veiling 
changes, the stellar brightness would change accordingly. 
However, there had been several hints that this is not al- 
ways the case. Gahm et al. (|2008l) demonstrated that in 
four cTTS with very pronounced emission line spectra 
and strong irregular light variations, the veiling variations 
clearly did not follow this pattern. Substantial changes in 
the veiling were accompanied by only minor changes in V 
and B, which in combination with an analysis of variations 
in equivalent widths of emission lines led the authors to con- 
clude that veiling is not solely related to an excess contin- 
uum emission. Furthermore, it was demonstrated that the 
photospheric absorption lines are weakened due to filling- 
in by narrow, emission line components. While the above 
shows that the veiling effect is caused by at least two dif- 
ferent processes, it is still unclear which mechanism causes 
the appearance of these narrow emission features. 

In order to investigate the origin of the narrow-lined 
chromospheric emission features, we performed a detailed 
investigation of the behavior of veiling in one of the most 
variable cTTS, namely DR Tauri. Our observations show 
this star to vary from a state of practically no veiling to 
levels of more than 10 times the continuum flux. At high 
veiling levels narrow emission features are clearly present in 
this object. In addition we find periodic variations in radial 
velocity of both absorption lines as well as of the narrow 
emission components, which, as we will show in this paper, 
provides indications on their origin. 

DR Tau is an unusually "active" cTTS, in the sense that 
the emission line spectrum is exceedingly rich and variable, 
and that the photometric variations are large and primar- 
ily related to physical changes in the plasma surrounding 
the star. There are many other well-known cTTS with sim- 
ilar properties, for which similar conclusions can be drawn 
regarding the nature of the veiling and the radial velocity 
fluctuations. We will illustrate this in the present paper by 
complementing our analysis with stars that show similar 
spectral line variations, namely RW Aur A, DI Cep, RU 
Lup, and S CrA SE. 

DR Tau was for a long period a relatively faint star 
until in 1976 it rather suddenly brightened by several mag- 
nitudes in the blue (Chavarria-K. I1979p . and the star is 
sometimes classified as an EXor star (see e.g. Lorenzetti et 
al. I2009p . However, since that time it has stayed at this 
bright level but with dramatic and fairly irregular light 
variations (Artemenko et al. I2010p . Tentative photomet- 
ric periods proposed by Richter et al. ^992) . Bouvier et al. 
(|1993p . and Bouvier et al. (|1995p differ by several days, and 
are different from a more recent determination by Percy et 
al. |2010| , who assign a period of 5.0 days. This is also in 
the same range of possible periods found from near-infrared 
(NIR) photometry by Kenyon et al. (jl994p . who relate the 
periodicity of DR Tau to a hot spot on a rotating star with 
a magnetic axis inclined to the stellar rotation axis. 

Pronounced spectral variability in DR Tau was recog- 
nized early (Bertout et al. 119771 Krautter & Bastian[19Mi 
Mundt ri984- Aiad et al. [T9841 Appenzeller et al. [T988|) . 
The strong emission lines on occasion show inverse P Cygni 
structures, which appear and disappear on time-scales as 
short as a few hours (Smith et al. 119971 119991) . Basri & 
Batalha (|1990p demonstrated that the veiling is large and 
variable, and subsequent investigations of relations between 



variations in emission lines and veiling led Guenther & 
Hessman (|1993p and Hessman & Guenther (|1997p to con- 
clude that in principle the observations were consistent with 
the concept of magnetospheric accretion, with gas flowing 
along dipole magnetic fields to a hot region on the vis- 
ible surface. However, the data are inconsistent with con- 
stant accretion in a simple magnetic dipole, because irregu- 
lar mass accretion and possible flare-line events distort the 
Balmer line profiles. In their long-term study of spectral 
variablity in DR Tau, Alencar et al. (|200ip were not able 
to recover any periodicity in the spectral changes. However, 
Hessman & Guenther (|1997p found evidence of a periodic- 
ity of 4.5-5 days in the equivalent widths of strong lines, 
and Johns & Basri (|1995p reported a 5.1 day period in com- 
ponents in Ha. Because these spectral periods are similar 
to some of the photometric periods presented above, one 
may suspect that both phenomena are related to rotational 
modulation, and that the period of rotation of DR Tau is 
in the range of 4.5-5 days. 

The complex emission line profiles and their dramatic 
changes can be better understood if the profiles are decom- 
posed into different components arising from different emit- 
ting volumes around the star. The accreting gas produces 
Balmer line emission, sometimes with inverse P Cygni ab- 
sorption components, and the broad metallic lines reflect 
the bulk motion of this in-falling gas (see e.g. Muzerolle et 
al. I1998al I1998bp . Normal P Cygni absorption components 
are seen as well, indicating outflow in the form of massive 
winds or jets. The P Cygni absorptions are particularly 
distinct in some far-ultraviolet lines, like the Mg ii doublet 
(Kravtsova & Lamzin HM^ Ardfla et al. [2002). The Hei 
lines, particularly the NIR line at 10830 A, display absorp- 
tion both from the in-falling and out-flowing gas as anal- 
ysed in a series of papers by Beristain et al. ( 2001 ) , Edwards 
etaL (HOnSl 2006), Kwan et al. (IMTfp and Fischer et al. 
(|2008p . In all, DR Tau stands out as a cTTS with an unusu- 
ally dynamic and very complex interplay between accreting 
and ejected gas. 

As in many cTTS, the He lines in DR Tau display a 
narrow emission component, which normally is attributed 
to post-shock gas at the footprints of the accretion flows 
(Lamzin et al ll996[ Beristain et al. I1998p . However, in 
DR Tau similar narrow components with a full-width half- 
maximum (FWHM) of ~ 20 km s^^ were also found in lines 
of Call, Fei and Fell (Batalha et al. 119961 Beristain et al. 
I1998p . In fact, Batalha et al. found a positive correlation 
between line strength and veiling and speculated that the 
narrow component is related to a heated stellar atmosphere 
close to the accretion zone - a "hot chromosphere" . 

Below we present our observations and reduction meth- 
ods in Section [S] The methods of analysis and results 
are given in Section [3] and are discussed in more detail 
Section |4l Section [5] summarizes our results. 

2. Observations and data reductions 

Observations of DR Tau were carried out during two runs in 
2007 and 2009/2010 at the Nordic Optical Telescope (NOT) 
with the fiber-fed echelle spectrograph FIES (Frandsen & 
Lindberg 2000). The data were processed with FIEStooQ. 
This software package is fine-tuned for reducing data ob- 
tained with the FIES spectrograph and performs the full set 

^ |http://www.not. iac.es/instruments/fies/fiestool/FIEStool.html 
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Table 1. DR Tau observations. See the main text for a 
more detailed description of the colmxms. 



Year 


HJD 245... 


N exp 


S/N ratio 


VFsiio 


avF 


2007 


4101.4056 


2 


65 


4.4 


0.5 




4103.5857 


2 


45 


1.5 


0.3 




4104.4278 


2 


55 


1.2 


0.3 




4105.6235 


2 


80 


7.5 


1.0 


2009 


5192.5623 


4 


110 


4.8 


0.5 




5193.4952 


3 


80 


7.0 


1.0 




5194.4984 


3 


60 


5.5 


0.5 




5195.4934 


4 


90 


3.8 


0.5 




5196.4748 


4 


70 


4.6 


1.0 


2010 


5197.5002 


4 


150 


7.0 


0.7 




5198.5539 


4 


90 


4.6 


0.5 




5199.5762 


2 


35 


5.0 


1.5 




5201.4676 


2 


60 


3.6 


0.5 



of echelle spectral reduction steps, including bias subtrac- 
tion, modelling of scattered light, extraction through two- 
dimensional modelling of the spectral orders, flat-fielding 
and wavelength calibration. 

The extracted spectra have a resolution of A/AA w 
47 000 and cover a wavelength range of 3650-7250 A. In 
Table [T] we summarize observation dates, the number of 
20-min exposures obtained in each night and signal-to-noise 
ratio per pixel in the continuum of nightly averaged spectra. 
The last two columns give the mean veiling factors at 6110- 
6160 A and their standard deviations (see Sect. 3.3.1). In 
addition, we obtained spectra of the K7V flare star HD 
28343 (BD -1-21° 652) for comparison purposes. 

This study focuses on DR Tau, but includes also sup- 
porting data from the two components of S CrA (Gahm 
et al. IMISl) . RULup (Stempels et al. 1^0071 . R W Au r A 
(Petrov et aL l200T|) . and DI Cep (Gameiro et aL 12006)) . 

3. Data analysis and results 

We start with a brief description of the spectra of DR Tau 
and the determination of the basic stellar parameters: tem- 
perature, gravity, radius, and rotation, and continue with 
an analysis of the veiling as a tracer of accretion, and radial 
velocities of narrow emission lines as indicators of chromo- 
spheric surface inhomogeneities. 

3.1. General spectral properties 

We present in Fig. [1] two typical segments of the DR Tau 
spectrum. Here, we averaged all 13 spectra to increase the 
signal-to-noise ratio. For comparison we also show a spec- 
trum of the template K7 star artificially veiled by a factor of 
four, so that the strength of the photospheric lines become 
similar to those of DR Tau. 

The main features in the spectrum of DR Tau are as 
described by Beristain et al. (P^998). The strongest emission 
is in the hydrogen Balmer lines and Can H and K. Weaker 
emission is seen in neutral and singly ionised permitted 
lines of Fei, Fell, Hei, Hen, Tin, Nai and Cai, and in 
forbidden [Oi], [Sn], and [Nn]. 

The weakest emission lines are narrow (FWHM ~ 20 
kms~^), while stronger emission lines of metals exhibit a 
two-component structure: a similar narrow component as 
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Fig. 1. Segments of the averaged spectrum of DR Tau (solid 
line) as compared to the K7 V template HD 28343 (thin 
line) veiled artificially by a factor of four. In this and other 
figures the wavelength scale is in the stellar rest frame. 
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Fig. 2. Emission line profiles with red-shifted absorption 
components in the spectrum of DR Tau on HJD 2455192. 
Left: H(5, right: Fen 5018.43 A (the narrow emission peak 
on the blue wing is Hei 5015.67 A). 



well as a broad component with FWHM ^ 100 kms""'^. The 
broad component gradually grows stronger as the emission 
lines themselves become stronger, as is evident from the 
panels shown in Fig. [1] The stronger Hei and Fen lines 
(for example Fen A 4924 and A 5018 A) show additional 
red-shifted absorption extending to -1-350 kms~^. These 
signatures of infall vary with time and are also present in 
the higher Balmer lines as shown in Fig. [51 The profiles of 
Ha and other Balmer lines are split by deep broad absorp- 
tion, centred at —100 kms"^, which indicates an outflow. 
Typical Balmer line profiles can be found in Alencar et al. 
(|200ip . where the Balmer emission was analysed in greater 
depth. 

Photospheric lines are detectable longward of - 4500 A 
in the averaged spectrum, while during some nights with 
higher veiling and / or lower signal-to-noise ratio many lines 
are not detectable. Nevertheless, in some spectral intervals 
the photospheric lines are always visible, which makes it 
possible to follow night-to-night variations of the veiling. 
In moments of high veiling some photospheric lines (e.g. 
Fei, multiplet 15) turn from absorption into pure emission. 
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3.2. Stellar parameters 

Spectral types assigned previously for DR Tau fall in a fairly 
narrow range, from K5 to MO. The variable and often high 
veiling of the photospheric absorption lines makes it diffi- 
cult to derive precise stellar parameters. However, we did 
obtain one spectrum of DR Tau in a relatively calm state 
(HJD = 2454104), which we will refer to as the low- veiling 
spectrum. 

This low-veiling spectrum was compared to a grid of 
synthetic spectra calculated for Tcff = 3500, 3750, 4000, 
4250, and 4500 K with log g values of 3.0, 3.5, and 4. 0. This 
grid was calculated using the code by Berdyugina (|199ip 
and Kurucz' models. Atomic line data were retrieved from 
the VALD database (Kupka et al. 1^000]) . The best fit was 
found for T^s = 4000 - 4250 K and log g = 3.5 - 4.0. In 
order to check the spectral type with molecular features, we 
used the spectrum of HD 28343 (spectral type K7V) along 
with MOV and K5V spectra collected earlier with FIES. 
The spectrum of DR Tau is best matched with the K7V 
spectrum, as is illustrated in Fig. [31 where the template 
spectrum is artificially veiled by a factor of 1.0. The TiO 
bandhead at 7054 A is also consistent with a spectral type 
of K7V for DR Tau, and we adopt T^s = 4100 ± 200 K and 
log g = 3.7 ± 0.5 for the photosphere of DR Tau. 

In Fig. SI we compare the low-veiling spectrum of DR 
Tau with a synthetic spectrum in the region of Cai A 
6124 A. The synthetic spectrum was calculated for Tcff = 
4100 K, log g = 3.7, Umic = 1.5kms~^, Umac = 1.5kms~^ 
and V sini = 5 kms~^, convolved with the instrumental 
profile of FIES and veiled by factor of 1.3. The profiles 
of the Vl, Til and Ban lines fit fairly well. A \ogg in the 
range 3.5-4.0 gives a reasonably good fit to the broad wings 
of the Cai line. The lines of Fel are partly filled in with 
chromospheric emission (see below). 

The projected rotational velocity v sini = 5 ± 1 kms^^ 
was found from the widths of many metal lines as fitted 
with the synthetic spectrum in the region 5580-6150 A. 
Note that v sini is somewhat less than the FWHM of the 
instrumental profile (w6.5 kms~^). Nevertheless, the high 
signal-to-noise ratio in our spectra enables us to detect this 
low value of the rotational broadening of photospheric lines. 

Assuming that at minimum brightness the photosphere 
is less affected by accretion, we can estimate the stellar lu- 
minosity and radius. From an extended set of photometric 
data from 1986-2003 (Grankin et aL lMITI) . Knin = 13'?^0± 
a°30 and the corresponding colour V -R = 1™40 ± ^"08, 
which gives Ay = 0™93 ± 0^"30. From the corresponding 
minimum B — V colours we find a faintest B magnitude 
that is similar to what DR Tau had in the first half of the 
20th century fChavarria-K. 11979)) . Hence, there is no indica- 
tion that the stellar parameters have changed following the 
general increase in brightness around 1976. The circumstel- 
lar component of the extinction is probably small because 
the star is seen close to pole-on (see also the discussion in 
Sect. EXll). With a distance to DR Tau of d = 140 pc (e.g. 
Bertout & Genova 2006) we obtain the following stellar 
parameters: L = 0.52 ± 0.15 and R = 1.46 ± 0.20 R©. 
Earlier published estimates of the radius range from 1.2 Rq 
(Bertout & Basri lT988|) to 2.7 R© (Johns-KruU & Gafford 

Regarding other programme stars, the primary and sec- 
ondary of S CrA have earlier been classified as G5 Ve and 
K5 Ve, respectively (Carmona et al. I2007p . However, our 



3 




7035 7040 7045 7050 7055 7060 7065 7070 7075 



Wavelength, A 

Fig. 3. Molecular TiO band at 7054 A. Upper: the low- 
veiling spectrum of DR Tau. Lower: the spectrum of 
HD 28343, (K7 V), artificially veiled by a factor of 1.0. 
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Fig. 4. Solid line: the low- veiling spectrum of DR Tau. Thin 
line: synthetic spectrum for Tcff — 4100 K, log g = 3.7 
and V sini = 5 kms^^, convolved with the instrumental 
profile and veiled by a factor of 1.3. Elements and lower 
level excitation potentials (eV) are indicated. 

long-slit UVES /VLT spectra yield very similar parameters 
for the two components, namely Toff — 4250 ± 150 K, \ogg 
— 4.0 ± 0.2 and v sini = 12 ± 1 kms~^. Furthermore, RW 
Aur A has a spectral type of K1-K4, v sini w 20 kms~^ 
(Petrov et al. IMJI]) : RU Lup has Tcff = 3950 ± 200 K, 
log (7 = 3.9 ± 0.3, and v sini = 9.0 ± 0.9 kms~^ (Stempels 
& Piskunov l2002p . and DI Cep has spectral type G8 and 
V sini = 20.0 ± 1.5 kms~i (Gameiro et al.HDOS)- 

3.3. Two sources of veiling 

As discussed early by e.g. Bertout (|1984p . the photospheric 
spectrum of a cTTS can be veiled by 1) an additional non- 
photospheric excess continuum and 2) line emission filling 
in photospheric absorption lines. 

3.3.1. Traditional veiling measurements 

Traditionally, the veiling factor (VF) is determined from the 
equivalent width (EW) of a photospheric line measured in 
the spectrum of DR Tau and in a template (e.g. synthetic) 
spectrum: 

VF - EWsynt/EWoR - 1. 

From measurements of numerous photospheric absorp- 
tions in the average spectrum of DR Tau no significant cor- 
relation was found either in veiling versus wavelength, or in 
veiling versus excitation potential of different lines, within 
the wavelength range 4500-6800 A and excitation potential 
range 0-5 eV. Although it is commonly accepted that veil- 
ing increases towards the blue, the wavelength dependence 
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Fig. 5. Ratio of EWs in the K7V template and DR Tau. 
Open circles: low-veiling spectrum, filled circles: average 
spectrum. The horizontal dashed line marks the case of 
zero veiling. 



is not always well defined because of the large scatter in 
the veiling factors derived from different spectral lines (e.g. 
Basri & Batalha. ri990|) . In Fig. [5] we show the ratio of EWs 
as measured in the spectra of the K7V template and DR 
Tau as a function of wavelength. Each symbol represents 
one spectral line. Note that the large scatter of points is 
caused not only by errors of measurements of EWs, but 
also by the dependence of veiling on line strength (see be- 
low). 

Not all absorption lines that are visible in the averaged 
spectrum can be detected in individual spectra. We there- 
fore carefully selected spectral regions where photospheric 
lines are clearly visible in all spectra of DR Tau, even in 
those with high veiling, and which never turn into emission: 
5580-5605 A (six lines of Cai), 6000-6040 A (four lines of 
Mm and Fei), 6IIO-6I6OA (five lines of Vi, Tii, Ban, 
Nai) and 6460-6480 A (t wo lines of Cai). In these regions 
we measured the corresponding values of the VF in each 
spectrum of DR Tau using the K7V template. Fig.[n]shows 
the night-to-night variations of the veiling. The VF at 6110- 
6I6OA and its standard deviation are also given in the last 
two columns in Table [TJ 

3.3.2. Distinguising the excess continuum from 
chromospheric line emission 

Unfortunately, the veiling measurements we perform above 
cannot discriminate between the effect of a veiling contin- 
uum and the effect of chromospheric line emission filling- 
in the photospheric lines. We therefore used the following 
three methods to separate and quantify these two effects. 

The first method is based on measurements of the for- 
bidden lines and is similar to the method used by Gahm et 
al. poos]) . Forbidden emission lines of [O i] at 5777.3, 6300.2 
and 6363.8 A and [Sii] at 6717.0 and 6731.3 A are usu- 
ally present in spectra of cTTS. These forbidden lines are 
formed in a large volume extending far away from the star, 
as compared to permitted emission lines, and can therefore 
be expected to be intrinsically stable on a time scale of days. 
Hence, the observed night-to-night variations of EWs of the 
forbidden lines are mainly caused by continuum variations. 

Fig. [7] shows that there is a strong correlation between 
EWs of photospheric absorption lines and forbidden lines of 
[O i] . The existence of this correlation indicates that both 
the photospheric and the [O i] lines are affected by the con- 
tinuous veiling component. However, the range of variabil- 
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Fig. 6. Night-to-night variations of the veiling factor in DR 
Tau. Filled circles: spectral region 6000-6040 A (1 ines of 
Mm and Fei). Open circles: spectral region 61 10-6160 A 
(lines of Vl, Til, Ban, Nai) 
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Fig. 7. Correlated variations of EWs of forbidden emission 
lines (sum of [Oi] A 5577, A 6300 and A 6363 A) and photo- 
spheric lines (sum of 17 lines within 5580-6500 A). Dashed 
line shows relation expected in case of only continuum vari- 
ations. 



ity for the photospheric lines is twice as large as for the 
[Oi] lines: the EWs of [Oi] lines change by a factor of 2.7, 
while the EWs of photospheric lines change by a factor of of 
5.3. This means that another source of veiling contributes 
as much as the continuum veiling component. 

The second method is based on the effect of differential 
veiling. Here, we measured from the average spectrum of 
DR Tau EWs of about 90 photospheric lines of different 
elements in the spectral region 5580-6720 A. We then mea- 
sured the same lines in the template spectrum. For each 
spectral line we calculated the ratio of EWs in DR Tau 
and in the template spectrum. Identical measurements were 
performed on the low-veiling spectrum of DR Tau. 

If the veiling were only caused by a featureless con- 
tinuous veiling component, all lines would be reduced by 
the same factor. However, the ratio of EWs depends on 
line strength: stronger lines are veiled more (see Fig. [5]). 
The effect as such has been recognized before (e.g. Basri 
& Batalha 1990 and references therein) but is here demon- 
strated in a graphical form. We obtained the same result 
when we used a synthetic template spectrum instead of the 
K7V data, and the effect of the differential veiling persists 
when temperature and gravity of the template vary within 
probable errors in Tcs (± 200 K) and log g (± 0.5). Note 
that the scatter of points in Fig. [S] is solely owing to un- 
certainties in the EW measurements from the spectra of 
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Fig. 8. Ratio of EWs of photospheric lines of DR Tau and 
the observed K7V template HD 28343 as a function of line 
strength. Filled circles - averaged spectrum of DR Tau. 
Open circles - low- veiling spectrum of DR Tau. Regression 
lines are shown for both sequences. The 99% confidence in- 
tervals are outlined by dashed curves. The scatter of points 
around the regression lines is caused by errors of measure- 
ments of EWs. The horizontal dashed line corresponds to 
zero veiling. 



DR Tau and the template. From the weakest lines (EW < 
0.05 A) in Fig.[8]we estimate VF ~ 0-0.5 for the low-veiling 
spectrum, and VF ~ 3-4 for the averaged spectrum. These 
numbers represent upper levels of the veiling continuum in 
the two spectra. The increasing excess of veiling towards 
stronger lines (EW > 0.1 A) should be attributed to chro- 
mospheric line emission. The VF values given in Table [1] 
and shown in Fig. [6] are derived from relatively strong pho- 
tospheric lines (EW « 0.1 A in the K7V template), because 
weaker lines are not always measurable in the individual 
spectra of DR Tau. 

The third method uses the colourimetric characteristic 
of a star with a hot spot (Petrov & Kozack|5DD7j) . According 
to standard magnetospheric accretion models the contin- 
uum veiling originates from the emission of a hot spot lo- 
cated below the accretion stream. The expected effective 
temperature of the hot spot is 6000-8000 K with a filling 
factor within 0.1-1 % of the stellar surface area (Calvet & 
Gullbring |1998p . It is believed that stochastic events of ac- 
cretion cause the irregular night-to-night variability of the 
brightness and colours of cTTS. 

The expected colours of a K7 star with a hot photo- 
spheric spot can be calculated as a function of spot size 
and temperature. These two parameters define the veiling 
factor VF, and hence the calculated colours. This means 
that, for instance, changes in B — V and V — R are di- 
rectly related to the degree of veiling. In order to analyse 
the colour changes in DR Tau we made use of the photo- 
metric catalogue by Grankin et al. (|2007|) . which includes 
about 20 years of observations of TTS. For DR Tau there 
are over 600 observations in 18 seasons from 1986 to 2003. 




0.0 0.5 1.0 1.5 a.O 0.0 0.5 1.0 1.5 2.0 
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Fig. 9. Two-colour diagrams with a grid of "star + spot" 
models, shown with the bow-shaped curves. The numbers 
along the curves represent the amount of veiling in each 
model. In the left panel, the cluster of dots are the observed 
colours of DR Tau. In the right panel, a correction for in- 
terstellar reddening Ay = l'."0 has been applied. The slope 
of the interstellar reddening is indicated with an arrow. 

The left panel of Fig. [9] shows colour-colour diagrams 
with the main sequence from AO to M6, as well as the ob- 
served colours of DR Tau. The two solid bow-shaped curves 
are sequences of the "star-|-spot" colours, corresponding to 
two different colour temperatures of the hot spot, approx- 
imately those of FO and AO stars. Along these curves, the 
spot radius varies, thus changing the filling factor and, cor- 
respondingly, the veiling factor VF. The loci of VF = 1, 2, 
and 4 are marked in the figure. The observed colours de- 
viate from the model grids, but corrections for interstellar 
reddening must be made. 

The right panel of Fig. [9] shows the same diagram with 
the colours of DR Tau corrected for an interstellar redden- 
ing Ay = 1™0. This correction brings the cluster of points 
to the sequence of the "star -I- spot" colours. The cluster of 
points is centred at the position of VF « 1.0, with a spread 
from VF — 0.5 to VF — 2. Some scatter towards bluer B- 
V colours (about 0™2) may be caused by contribution of 
emission lines to the B pass band. 

Apparently, the hot region is permanently present on 
the visible surface of the star, but does not contribute to 
the veiling by more than VF « 2.0. This corresponds to a 
brightness variability in the V band within « 1™.2. Because 
the observed veiling of DR Tau reaches much larger val- 
ues, typically 3 to 10, the rest of the veiling effect must be 
caused by the line emission fiUing-in the photospheric lines. 
Generally, if the effect of the line emission reduces the pho- 
tospheric line depth by a factor of VFune and, in addition, 
the non-photospheric continuum reduces the line depth by 
a factor of VFcont, the resulting veiling VF is defined by 
the relation 

(VF + 1) = (VFiinc + 1) • (VFeont + 1). 

This is a steep relation. For instance, with VFiinc — 
VFcont = 2.0, the resulting veiling is VF = 8. This explains 
the extremely high veiling (VF > 10) observed sometimes 
in spectra of cTTS. It also means that the accretion rate de- 
rived on the assumption that the observed veiling is caused 
by an additional continuum may be overestimated by a con- 
siderable factor, as was also concluded earlier by Gahm et 
al. (2008) for other cTTS. 

Note that the temperature used in this method is a 
colour temperature in the spectral range covering the BVR 
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Fig. 10. Anti-phase RV-variations of photosphcric (filled 
circles) and clironiospheric (open circles) lines in DR Tau. 



pass bands, and may differ from the effective temperature. 
The continuum radiation from a hot spot on a cTTS is 
composed of an optically thick continuum from the heated 
photosphere and an optically thin Paschen continuum from 
the pre- and post-shocked gas (GuUbring et al. I2000p . 

We conclude that all three methods described above 
show that the veiling continuum in DR Tau is relatively 
small compared to the values derived directly from mea- 
surements of EWs of the photosphcric lines, indicating that 
the excess of the observed veiling is caused by the chromo- 
spheric line emission. 
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Fig. 11. Anti-phase RV-variations of photosphcric (filled 
circles) and chromospheric (open circles) lines in four cTTS. 
Narrow emission lines of Fe i were measured in RW Aur and 
S CrA SE, and of He i in RU Lup and DI Cep. Note the 
positive shift in velocity of the He i emission. 



3.4. Anti-phase radial velocity variations 

Stellar radial velocity (RV), as measured from positions of 
photosphcric lines, can vary on time scales of days either be- 
cause of orbital motions in a close binary system or because 
of surface brightness inhomogeneities (cool/hot spots). In 
the latter case variations in RV must be within the observed 
V sinj value. 

Stellar RVs of DR Tau were measured by cross- 
correlation with the spectrum of the template HD 28343 
in the region 5560-6670 A, with the emission lines masked 
out. The radial velocity of the template was determined by 
cross-correlation with the synthetic spectrum {Tes = 4000 
K, log g = 4.5) resulting in RV = -35.20 ± 0.04 kms"^ 
The results of the RV measurements are given in Table [5J 

We then selected a group of unblended emission lines 
of the Fei multiplet 15, within 5371-5456 A, as representa- 
tive of the chromospheric emission (most lines are present 
in Fig. [1]). The profile of each line was decomposed into a 
broad and narrow component, and the position of the nar- 
row component was determined by a Gaussian fit. The av- 
erage RVs of the Fe i emission lines are also given in Table [2] 
The two blanks in the column of RVcm correspond to dates 
when the star was at a state of low veiling and the Fe i emis- 
sion was too weak to be used for RV measurements. The 
RVs show night-to-night variations with a semi-amplitude 
of about 1 kms~^ for both the photospheric and the emis- 
sion lines, which is well within the v sin i value. The main 
result is ttiat tlie RV variations of the photospheric lines and 
of the narrow component of the Fe i emissions are in anti- 
phase (see Fig. [T0|) . Similar variations were first reported 
for RW Aur by Petrov et al. 

The broad components show periodic shifts in central 
velocity with time. However, the measurement accuracy of 
the position of a broad component using a Gaussian fit is 



much lower than for the narrow components. The RV am- 
plitudes of broad components were measured by means of 
cross-correlation of each spectrum with the average spec- 
trum of DR Tau with the narrow components masked out. 
The amplitudes of the broad components turned out to be 
larger than those of the narrow components, which is con- 
sistent with an origin of the broad components well above 
the stellar surface. 

Table [3] summarises measurements of RVs in different 
groups of lines, which are arranged in order of increasing 
RV amplitude. The mean stellar velocity is -|-23 kms"""^. 
Note that the mean RV is more positive for the "high- 
temperature" lines, up to +36 kms^^ for the narrow Heii 
component. 

Our search for radial velocity changes was then ex- 
tended to the other cTTS included in the present study. As 
shown in Fig. llll all programme stars show the same type of 
anti-phase variations as observed for DR Tau. Furthermore, 
narrow emission lines from ions like Hei formed at higher 
temperatures are also red-shifted by several kms^^ with 
respect to the stellar velocity. Therefore these phenomena, 
indicating infall to regions at the stellar surface that are off- 
set relative to the rotational axis, appear to be a universal 
property of active cTTS. 

4. Discussion 

4.1. Modelling the RV variations 

Cross-correlation of DR Tau spectra with the synthetic 
spectrum in the region 6000-6040 A reveals that the bis- 
sector of cross-correlation function (as described in Dall et 
al. 2006) correlates with the stellar radial velocity. This 
excludes the case of a binary motion and suggests that RV 



8 P. P. Petrov et al.: Accretion-powered chromospheres in classical T Tauri stars 

Table 2. Radial velocities of photospheric and chromospheric lines in DR Tau 
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Table 3. Mean RV and range of RV variations 



Line 


Mean RV 


Semi-amplitude 




(kms-i) 


(kms-^) 


Photospheric lines 


23.0 ±0.2 


1.1 ±0.2 


Narrow component: Pel 


22.5 ± 0.2 


1.1 ± 0.2 


Fell 


23.4 ± 0.3 


1.5 ± 0.2 


He I 


26.9 ± 0.9 


1.1 ± 0.3 


Hen 


36.1 ± 0.5 


2.5 ± 0.7 


Broad component: Pel, Fell 


24.3 ± 1.0 


5.0 ±1.5 


He I 


22.9 ± 1.0 


8 ±2 



variations are caused by distortions in the photospheric line 
profiles. The same effect was found earlier for RU Lup by 
Stempels et al. ([MI71) . 

The anti-phase RV variations of the photospheric and 
emission lines suggest that the line emission is confined to a 
certain area on the stellar surface, which is not aligned with 
the axis of rotation, resulting in periodic shifts in RV as the 
star rotates. The velocity curves in Fig. 1101 have minima 
and maxima separated by about 4.5 days, which is similar 
to the photometric and spectroscopic periods (4.5-5.0 days) 
reported above, and which have been linked to rotational 
modulation. 

The effect of an "emission spot" resulting in the anti- 
phase RV variations can be illustrated with a simple model. 
We generated synthetic spectra of a single rotating star plus 
an area emitting the narrow emission lines under the as- 
sumption that this area has the form of a circular spot. The 
input parameters for the model are: projected rotational ve- 
locity {v sin i), inclination (i) of the stellar rotational axis to 
the line of sight, spot latitude (0), spot radius (rgpot), local 
spectra of the photosphere and the chromospheric spot, and 
the limb darkening coefficient (/x). The spectrum was inte- 
grated over the visible stellar surface at different phases of 
stellar rotation. We then determined the radial velocities of 
the photospheric lines RVphot and the emission lines RVcm 
as a function of phase. 

In this model most of the input parameters are re- 
stricted by observations. From v sini = 5 ± 1 kms~^, ro- 
tational period 4.5-5.0 days and R* = 1.46 Rq, we obtain 
the inclination angle i = 20° ±4°, i. e. the star is seen 
close to pole-on. The inclination is not very sensitive to the 



Table 4. Stellar and chromospheric spot model parameters 



Stellar parameters 





log 9 


L, 




V sin i 


i 


(K) 




(L«) 


(Ro) 


(kms^i) 


(deg) 


4100 


3.7 


0.52 


1.46 


5.0 


20 



Chromospheric spot parameters 



^spot 






(deg) 


(deg) 




14 


76 





adopted Ay, because the slopes of lines of equal radii in 
the V versus {V — R) diagram and reddening are the same. 

The spot latitude ((/>) on the stellar surface can be esti- 
mated from the amplitude ARV of the radial velocity varia- 
tions of the Fei emission lines: ARV= wsini • cos(0). With 
vsini^ 5.0 ±1.0 kms'^ and ARV:= 1.15 ±0.20 kms'^ we 
get (j) = 76° ±4°. 

A synthetic spectrum with Tcs — 4100 K and log 5 = 3.7 
was taken as the local spectrum of the undisturbed photo- 
sphere surrounding the spot. Within the spot area the local 
spectrum was assumed to be a sum of the photospheric 
spectrum and a line emision spectrum with similar line 
widths and line ratios. The overall emission line strength 
was a free parameter. The limb darkening coefficient was 
taken as ^ = 0.8 (Al-Naimi ll978p . Table S] summarises the 
model parameters. 
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Fig. 12. Upper panels: location of the chroniospheric spot 
and the local spectra corresponding to the undisturbed pho- 
tosphere (absorption) and to the spot (emission). Lower 
panels: left - profiles of the Cai A 6039 A line, correspond- 
ing to phases 0.0 spot facing the observer, solid line), 0.22 
and 0.78 (maximal blue- and red-shifted absorption, dotted 
lines). Dashed profile: a normal photosphere without spot; 
right - RV- variations of absorption lines (open circles) and 
emission lines (filled circle). 

As shown in Fig. [T2] the model can reproduce the ob- 
served anti-phase RV- variations very well. The emission line 
strengths were tuned to obtain the best fit to the observed 
RV amplitudes. Equal RV amplitudes of emission and ab- 
sorption lines occur when the latter is about halfway filled 
with emission (i.e. VFunc ~ !)• A spot radius of rgpot — 14° 
(filling factor = 0.015) adopted in this model is about the 
same as the typical radius of the accretion stream in the 
3-D simulations of magnetospheric accretion by Romanova 
et al. (|2003p . The model also reproduces the redshifted ve- 
locity offset of the He lines, provided that this emission 
is formed in infalling gas, with a velocity vector directed 
towards the star. 

In some other cTTS, the Zeeman-Doppler imaging tech- 
nique has been used to reconstruct the large-scale magnetic 
topology and to outline Ca ii line emission from the shocked 
gas near the stellar surface: V2129 Oph (Donati et al. 120071 
and [MH), BP Ta u (Do nati et al. ,2008). CV Cha a nd CR 
Cha (Hussain et al HUU^ . and AA Tau (Donati et al. [2(mi| . 
While the recovered structures are sometimes complex, it 
is not uncommon to see hot spots, as indicated by the en- 
hanced Can emission. The adopted size and derived loca- 
tion relative to the rotational axis for our Fe i emitting area 
in DR Tau resembles the spots of Can emission found in 
V2129 Oph and BP Tau. 

It is worthwhile to note that the distortions of a photo- 
spheric line profile caused by an emission core (Fig. I12p are 
very similar to those caused by a dark spot. The two effects 
can be distinguished with the help of photometry made si- 
multaneously with spectropolarimetry. Although we can- 
not discard a contribution from cool spots in the observed 
RV-variations of photospheric lines in DR Tau, both the 
differential veiling and the anti-phase RV-variations indi- 
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Fig. 13. Correlated variations of RV of the narrow emission 
of Fei and Hen. 

cate that these are primarily caused by line emission. The 
presence of narrow emission lines inside photospheric lines 
is particularly evident in S CrA SE (Gahm et al. 12008^ . 

4.2. Accretion powered emission 

A direct signature of infall are the inverse P Cyg (IPC) ab- 
sorption components. In our series of spectra of DR Tau, 
IPC are present in the strongest Fen lines (multiplet 42), 
in the Nai D lines, in the higher Balmer lines and, less evi- 
dent, in the He I line at 5876 A (Fig. (2). In the low- veiling 
spectrum most emission lines have faded out, and only the 
strongest emission lines are present. On this occasion the 
IPC profiles have also disappeared. During this episode of 
relatively weak accretion both the veiling continuum and 
the chromospheric emission are considerably reduced. The 
upper limit to the veiling continuum that day was only 
0.5, as is evident from Fig. [Sj The area of enhanced chro- 
mospheric emission in DR Tau is permanently present, but 
varies in strength, and is probably coupled to a variable 
amount of mass accretion. 

We can roughly estimate the expected accretion lumi- 
nosity for the low-veiling spectrum using the relation be- 
tween Ha emission line flux and accretion luminosity by 
Fang et al. (12009 1. Our spectra are not flux-calibrated, 
but for the low-veiling spectrum we can assume that DR 
Tau was close to its minimum brightness (see Sect. 3.2). 
In the low-veiling spectrum, EW{Ha) — 83.5 A. For a 
K7V star with = 0.5 Lq, this EW corresponds to 
log{LHa/ Lq) = —2.53. Using the relation by Fang et 
al. (,2009) . we obtain log{Laccr/ Lq) = -0.89 ± 0.25, or 
Laccr/Li, ss 0.26, which is compatible with an observed 
veiling < 0.5 in the V band. 

To explore further whether the narrow Fe i emission is 
related to the hot accretion spot, we followed radial velocity 
variations of the He n A 4686 A emission. In standard dipole 
magnetospheric accretion models He is ionised by the X-ray 
radiation from the shock. Fig. [13] shows that the RV of the 
narrow Fe i emission lines vary in phase with the narrow 
He n emission line. This strongly suggests that the area of 
chromospheric emission is located close to the hot region at 
the footprint of the accretion stream (see also discussions 
by Batalha et al. '1996'and Beristain et al. ll998|) . Moreover, 
the Hen line is clearly red-shifted relative to other nar- 
row lines. As discussed by e.g. Lamzin (|1998p . these "high- 
temperature" lines are formed in the post-shock region. 
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Our observations thus demonstrate that there is a 
close relation between the appearance of enhanced chromo- 
spheric emission and accretion in classical T Tauri stars. 
The narrow emission lines of e.g. Fe i appears over an area 
at the stellar surface that is associated with the footprints 
of the magnetic funnels, and that is offset from the pole of 
rotation. 

In standard dipole magnetospheric models the accretion 
spot may take the form of a ring, which becomes distorted 
when the magnetic and rotational axes are not aligned. 
However, recent models assume much more complex mag- 
netic field configurations (e.g. Mohanty & Shu '2008). The 
3-D simulations of magnetically channelled accretion by e.g. 
Romanova et al. (|2004l [20TT|) and Long et al. pOTT|) reflect 
the complex magnetic structures found from Doppler imag- 
ing of TTS. 

In the general concept of magnetically channeled accre- 
tion, the accretion spot on a cTTS determines the observed 
brightness of the star, the veiling of photospheric lines and 
the strength of certain emission lines, and all three param- 
eters must vary in correlation. However, observations do 
not always support this prediction (see e.g. Petrov et al., 
I2001|l . In order to resolve this discrepancy, we suggest that 
the area of enhanced chromospheric emission is much more 
extended than the size of the hot accretion spot. In the 
models by Romanova et al., less dense and non- uniform 
gas surrounds the main accretion stream and falls onto a 
much larger area than what is covered by the main stream, 
up to 16% of the stellar surface. A possibility is that our 
observed enhanced chromospheric emission is excited by 
such a milder accretion flow and therefore distributed over 
a similar area surrounding the smaller hot spot. One can 
assume that the accompanying shocks and related phenom- 
ena are less pronounced than at the main footprints, and 
that little or no continuous excess emission is released. In 
this scenario, a tight correlation between the brightness, 
veiling and emission line strength is not expected. 

X-ray observations of cTTS show that besides very pow- 
erful coronal emission there is as a rule one component of 
soft X-rays that can be associated with the hot spot (see 
e.g. Giidel & Telleschi [20071 and Argiroffi et al [2U09l and 
references therein). Recently, Brickhouse et al. (|2010p pre- 
sented X-ray spectroscopy of the cTTS TW Hya. A low- 
temperature component (T ~ 2.5 MK) could be identified 
with the accretion shock, as depicted in the standard ac- 
cretion models, and a high-temperature component was as- 
cribed to coronal emission at T ^ 10 MK. In addition, 
the authors discuss a third component, represented by low- 
density lines of e.g. O vii formed at T ^ 1.75 MK, which fills 
a great volume outside the post-shock region. Brickhouse 
et al. speculate that violent mass outflows and shocks may 
propagate from the base of the accretion shock (see also 
Orlando et al. I2010p , and that the photospheric surround- 
ings are heated, which supplies ionised gas to accretion-fed 
magnetic structures, such as magnetic loops. In addition, 
one can expect that the surface layers surrounding the hot 
spot are excited by UV- light directly from the hot spot. 

Regardless of the dominating cause of the enhanced 
chromospheric emission one would expect that the atmo- 
spheric structure adjacent to the hot spot is affected. In 
some early attempts to explain line emission and veiling in 
cTTS, Herbig (fTOTOl) and Cram f|1979. 1980) introduced the 
concept of a " deep chromosphere" , where the temperature 
minimum occurs at deeper layers in the atmosphere, thus 



resulting in the appearence of stronger emisson spectrum. 
We speculate that our area of enhanced emission surround- 
ing the hot spot could be similar to a deep chromosphere 
with a temperature minimum much deeper in the atmo- 
sphere than in a normal dwarf, but still above the photo- 
spheric continuum. In any case, the chromospheric emission 
in cTTS appears to be related not only to solar-like mag- 
netic activity, but is to a considerable extent powered by 
accretion processes. 

5. Conclusions 

We have studied periodic changes in photospheric lines and 
narrow components of emission lines in DR Tauri and other 
classical T Tauri stars, and investigated the nature of the 
dilution of photospheric lines, the so-called veiling. In par- 
ticular we note: 

— In the course of our observations the amount of veiling 
in DR Tau varies from practically nothing to factors 
of more than 10 times the stellar continuum intensity. 
We have developed three methods to distinguish the 
contribution from line emission and continuous excess 
emission to the veiling. We conclude that narrow emis- 
sion lines fill in the photospheric absorption lines, and 
that this effect is more pronounced when the veiling 
is strong. Therefore, estimates of mass accretion rates 
based on the assumption of a continuous veiling will be 
overestimated for cTTS with strong emission line spec- 
tra. 

— Radial velocities of photospheric and narrow chromo- 
spheric lines vary in anti-phase in all cTTS studied here. 
The mean velocity of the narrow lines are at rest rela- 
tive to the stellar radial velocity. The velocity ampli- 
tudes are small. We have modelled these variations and 
infer that the emission lines are formed in an area at 
the stellar surface that is offset from the rotational axis. 
These properties appears to be a common phenomenon 
in cTTS. 

— The narrow emission line of He ii A 4686 A in DR Tau 
varies in phase with the narrow chromospheric compo- 
nents, but is systematically red-shifted by ^ 10 kms~^, 
which is consistent with the origin of this line in de- 
celerated post-shock gas. We conclude that the area of 
enhanced chromospheric emission is directly related to 
the hot spot area, and that chromospheric emission is in- 
duced by accretion. In accordance with this, we demon- 
strate that a sudden disappearance of inverse P Cygni 
profles in strong lines coincided with a fading of both 
the veiling and the strength of chromospheric lines. 

We suggest that the area of enhanced chromospheric 
emission comes from an extended region at the stellar sur- 
face surrounding the hot spot. The physical mechanism that 
triggers the emission can be related to current ideas that 
differ from the standard model in that accretion is more 
wide-spread, or emanates from injection of mass and energy 
from the shocked region under the main accretion stream. 
As a result the photosphere surrounding the hot spot is 
heated, which leads to a modified atmospheric structure, 
which in turn produces the enhanced emission in narrow 
emission lines responsible for the extra component in the 
veiling through line-filling of photospheric absorption lines. 
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